clear
% 定义测试参数
nx0 = 64; % 示例中x方向的内部网格点数
ny0 = 64; % 示例中y方向的内部网格点数
alx=1;
aly = 5;  % y方向的长度，用于计算频率相关的部分
%dx=alx/(nx0+1);
% 调用函数生成矩阵A3
A3 = generate_poisson_matrix(nx0, ny0, alx,  aly); % 注意：alx在此处被替换为dx，根据您的具体应用场景调整
